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Abstract 

Pulsars, Gamma Ray Bursts and gamma ray emitting Active Galac- 
tic Nuclei represent the most luminous compact gamma ray sources in 
the Universe. Based on their apparent luminosities and non-thermal 
spectra, a direct association of the gamma ray sources with objects of 
the size of a neutron star or an accreting black hole appears question- 
able at a first glance: Compact radiation fields are optically thick to 
pair creation. The paradox can be resolved, if the gamma ray emit- 
ting plasma is in relativistic bulk motion with respect to the observer. 
The relativistically expanding plasma dissipates a significant fraction 
of its bulk kinetic energy into gamma rays, presumably through par- 
ticle acceleration at coUisionless shocks. Details of the conversion of 
bulk kinetic energy, carried by baryons, into radiation, emitted by 
leptons, are highly controversial. Results are important for the prob- 
lem of the origin of the highest energy cosmic rays and the compo- 
sition of the plasma emerging from the ergosphere of compact ob- 
jects. Measurements of ultra-high energy neutrinos can be used to dis- 
criminate between baryon-induced electromagnetic cascade or lepton- 
induced inverse- Compton models. 

Gamma rays from remote quasars at high redshifts may be used to 
probe the density of diffuse intergalactic infrared-to-ultraviolet pho- 
tons from the era of galaxy formation. An exponential turnover in 
the gamma ray spectra marks the onset of pair creation by photon- 
photon interactions between the quasar gamma rays and diffuse low 
energy photons along the line of sight. If the background density can be 
inferred from direct observations, one could determine the Hubble con- 
stant at unprecedented distances where distortions of the cosmological 
expansion due to mass fiows accelerated by local mass concentrations 
can be neglected. 



submitted to Rev.Mod.Astr. 

"Thus we conclude to attribute the steady-state to homology, kinetics however to 
anomaly. Anisotropy, on the other hand, is the reason for the anomalous nature." 



1 Introduction 



Ludwig Biermann's (1951) discovery of the interaction between the corpus- 
cular radiation from the sun and the ions in the tails of comets has forcefully 
demonstrated the importance of dynamical processes in the solar system. 
Satellite experiments demonstrate that the interplanetary medium experi- 
ences a wealth of phenomena leading to particle acceleration and non-thermal 
emission. Since the solar wind is rather slow, these non-thermal processes 
are only of local importance. However, the stellar wind of rapidly rotating 
stars and the ejecta of supernovae reach much higher velocities and, con- 
sequently, their non-thermal emission is much stronger. Radio synchrotron 
emission from the expanding shell of a supernova indicates the acceleration 
of particles to relativistic energies at the blast wave. Supernova remnants are 
the most likely origin of the cosmic ray particles below ~ 10^^ eV (Axford 
1981) or, possibly, to even higher energies (Volk & Biermann 1988). The 
presence of a frozen magnetic field in the highly conducting plasma intro- 
duces anisotropics in the phase-space distribution of particles driving them 
away from thermodynamic equilibrium. Energetic particles in the low den- 
sity plasma cannot maintain detailed balance by rare Coulomb collisions. The 
gyrating particles are forced to phase-space diffusion through a resonant cou- 
pling with the collective degrees of freedom of the magnetic field fiuctuations 
which are correlated by turbulent cascades. Today we know that stellar winds 
are a common phenomenon, and often intimately linked with the presence of 
rotation and magnetic fields. Parker (1958) has realized that the lever arm 
provided by the corotation of the stellar wind below the Alfven point along 
each streamline acts as a breaking torque extracting rotational energy and 
angular momentum from the star. It is therefore not surprising that rapidly 
rotating magnetized neutron stars and accretion disks around black holes 
produce powerful magnetized winds. Rotating disk winds can self-coUimate 
attaining a bipolar topology as found in radio jets emerging from many AGNs 
and some galactic black hole candidates (Camenzind 1990). 

In this Lecture, I will introduce the audience to various phenomena of 
gamma ray astrophysics and their possible explanation as the consequence 
of particle acceleration in the plasma fiows emerging from compact objects, 
viz. Active Galactic Nuclei (AGNs), Gamma Ray Bursts (GRBs) and pul- 
sars. Emphasis of the Lecture is on AGN jets. Experimental progress in 
this field of astrophysics has been very rapid in recent years. Astronomical 
objects can now be observed from radio waves to TeV gamma rays using a 
diversity of telescopes, from radio antennas to Cerenkov mirrors and scintilla- 
tor arrays providing sensitivity for 20 orders of magnitude in photon energy. 
More than a thousand GRBs (BATSE 3B catalog) and a surprisingly large 
number of AGNs (~ 50) have been detected by the Compton Gamma Ray 
Observatory (Eichtel et al. 1994, Thompson et al. 1995). Using ground-based 
Cerenkov telescopes, steady TeV (10^^ eV) emission has been discovered from 
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two pulsars, the Crab nebula and PSR 1706-044 (Weekes 1994), and from two 
AGNs, Mrk401 (Punch et al. 1992) and Mrk501 (Quinn et al. 1996). The 
gamma ray emitting AGNs belong to the class of blazars characterized by 
drastic variability at all wavelengths, a high degree of linear polarization and 
a broad-band continuum spectrum. Morphologically, the sources show com- 
pact radio cores and often jets. A common spectral characteristic of blazars, 
GRBs and the pulsed phase of pulsars (especially old pulsars) is that the 
energy flux is typically dominated by gamma rays (Montigny et al. 1995). 
Models of the radiation processes responsible for the ultraluminous gamma 
ray emission allow for a diagnosis of the plasma conditions in the vicinity of 
the compact objects powering the sources and along the line of sight. 

Gamma ray astrophysics brings together astronomy and physics on sev- 
eral important topics. Stable supersymmetric particles have escaped their 
detection in accelerator experiments because of their large masses, but they 
could be detected through gravitational interaction on astronomical scales. 
Their existence could solve the long-standing mystery of the missing mass 
inferred from flat galactic rotation curves or galaxy cluster dynamics. Direct 
evidence could be obtained by measuring high-energy neutrinos from dark 
matter particle-antiparticle annihilation in astronomical bodies where dark 
matter accumulates over the Hubble time (Kamionkowski 1992, Gaisser et al. 
1994). Astrophysical gamma rays indicate the possible existence of a diffuse 
flux of high-energy neutrinos which would constitute a conservative back- 
ground for searches of such 'new physics'. Astronomy provides many more 
rather challenging phenomena than just dark matter. There is mounting ev- 
idence for black holes, with the best evidence to date for the supermassive 
black hole in NGC4258 inferred from the Keplerian rotation of water-maser 
emitting clouds (Miyoshi et al. 1995). It is an unsolved problem in physics 
whether the existence of the event horizon, which develops during a gravita- 
tional coUaps, implies the existence of a space-time singularity in its interior. 
By Penrose's theorem, a singularity must exist if the general theory of rel- 
ativity is true. Only the yet undiscovered laws of quantum-gravity could 
forbid the formation of a singularity which we render unphysical based on 
our scientific belief. In a sense, the existence of black holes is telling us that a 
theory of quantum-gravity must exist. Astronomical observations can bring 
us very close to the event horizon in accreting black hole systems where mat- 
ter is dissipatively sinking into the potential well of the black hole. As Ludwig 
Biermann's corpuscular radiation from the sun, these black hole systems eject 
high speed particles along their angular momentum axes. The nature of these 
particles is still poorly understood, and in this Lecture I will present diagnos- 
tic strategies based on gamma ray and neutrino astronomy. Another exciting 
problem comes from cosmic ray astrophysics where trans- Greisen particles 
with energies exceeding 5 X 10^^ eV have been discovered recently. Based on 
the shower properties, the particles are most probably protons (Dedenko et 
al. 1995). Due to photo-production of pions in the 3K-background radiation 
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field, cosmic ray protons above the Greisen-energy can not propagate more 
than ~100 Mpc (Greisen 1966, Zatsepin & Kuz'min 1966). Where do the 
protons come from? They could be related to topological defects predicted by 
inflationary cosmology due to a phase transition in the early Universe (Sigl 
et al. 1994). However, new physics is not necessarily required by these exotic 
events, if protons are shock-accelerated in the jets of radio galaxies (Biermann 
& Strittmatter 1987). 

The plan of the Lecture is as follows: In Sect. 2 I will confront the audience 
with the striking paradox that the observed gamma ray point sources are too 
compact to emit gamma rays by introducing the important process of photon- 
photon pair creation. It is shown that relativistic bulk motion resolves this 
paradox. In Sect. 3 I will outline a scenario of accreting supermassive black 
holes driving a relativistic jet as the central powerhouse in AGNs. Sect. 4 will 
zoom-in on the proton-induced cascade mechanism for the gamma ray emis- 
sion in AGNs. In addition, flux predictions for ultrahigh energy neutrinos are 
presented and compared with experimental sensitivities in Sect. 5. A natural 
source for neutrinos with energies exceeding the electroweak phase transition 
energy would unlock the door to a new era of particle astrophysics. Finally, 
the notion of a gamma ray horizon due to intergalactic absorption of gamma 
rays is introduced in Sect. 6. Measurements of gamma ray absorption may 
provide important clues on the cosmic distance scale and on galaxy forma- 
tion. In addition, a strategy is developed to solve the halo vs. cosmological 
origin problem for GRBs. 



2 The Paradox 

Let us briefly look at the gamma ray luminosities and sizes of pulsars, GRBs 
and AGNs in Tab.l. The luminosities given in the first column correspond 
to the isotropic luminosity 

L = i.D^ E^dE (1) 

where D denotes the distance (for GRBs, see discussion after Eq.(13)), E the 
observed photon energy in the gamma ray band between -Emm and -Emax for 
a photon number distribution dN/dE = No[E / Eo)~^"~^^\ The sizes in the 
second column correspond to the neutron star radius Tq ~ ^hc/Gm'^h /nipC ~ 

3 km or the Schwarzschild radius Tq = rs = 2GM/c^ = 3 X lO^^mio cm of 
an Eddington radiating black hole of mass M = lO^^mioAf©, respectively. 
Column three shows the upper limit on the source size Tyar < cAt due to 
the requirement that the light travel time across the source is less than the 
observed variability time scale"^ At. The fourth column denotes the optical 

^Incoherent emission is tacitly assumed here. If coherent emission is important, the 
causality constraint would be invalidated. However, the size of phase-space is cx 
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depth to pair creation at a photon energy of = 2meC^ ~ 1 MeV (Gould & 
Shreder 1966), i.e. the process 

7 + 7^e+ + e- (2) 

Pair creation can only occur, if the two interacting photons provide more 
than twice the electron rest mass in their center- of- momentum frame, that is 
if the target photon energy e exceeds the threshold energy eth 

^eth(l + zf{l - cos = 2{m,^ f (3) 

in a collision at an angle 6^- The {l-\-zY factor takes into account cosmological 
or gravitational redshift. The interaction cross section is given by 



16 ^ ' 



2/3(/32-2) + (3-/3^)ln 



(4) 



where /3 = y 1 — Ijv? with = e/eth, and where ctt denotes the Thomson 
cross section. For an average collision in an isotropic target radiation field 
(cos = 0), the cross section Eq.(4) peaks at ~ (Jt/S when E ~ 2(meC^)^/e. 
The optical depth r-^-^ for gamma rays of energy E = 2meC^ is then given by 



T^^^o ^ n^{> meC^)-aTr (5) 

which is of sufficient accuracy for the qualitative comparison in Tab.l. An 
inverse power law distribution of target photons n^[> e) = dn oc e~" leads 
to the energy-dependent optical depth 

T^^{E) = T^^,o{E/2m,c^)" , (6) 

since the threshold condition Eq.(3) requires that the target photon energy 
obeys eth ^/E^ so that n^{> eth) oc ^th ^ E" . 

From Tab.l it is evident that the optical depth for E > 2meC^ exceeds 
unity implying the paradoxical result that the objects cannot emit non-thermal 
gamma rays above MeV. 

Pulsars: For the estimate of the optical depth of the Crab pulsar mag- 
netosphere, I have used the radius Tq = 10 km instead of cAt which is more 
appropriate for the emission from the pulsar wind. The pulsar wind has a 
very low gamma ray optical depth and there is no paradox for the unpulsed 
emission from the nebula. However, for the pulsed component the paradox 
becomes striking considering the energy dependence in Fq.(6). Pulsed gamma 
rays have been observed, e.g. by FGRFT at photon energies of 10 GeV yield- 
ing T-y-y ~ 100 (ctx = 0.5). An obvious way out of the problem is to assume 



making it very unlikely that coherent gamma ray processes occur in nature. 
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Table 1: Apparent pair creation optical depths for various gamma ray sources 



Type of Source 


L [erg/s] 


To [cm] 


cAt [cm] 


'^-77,0 


To 


Pulsar (Crab") 


10^6 




(10^^) 


1 


1 


AGN (3C279) 


^048 






100 


3 


GRB 




10^ (?) 


10« 


10^3 


100 



"pulsed component 

that the emission is not isotropic. In the case of a pulsar we know this inde- 
pendently by the very fact that we observe pulsed emission from a rotating 
beam. If the solid angle subtended by the beam is small, i.e. <^ 47r, the 
power emitted is only (O/47r)10"^^ erg s~^. Additionally, the photon-photon 
scattering angles are constrained to a narrow band increasing the threshold 
energy Eq.(3) and decreasing the pair production cross section Eq.(4). Thus, 
the paradox is removed. The physics responsible for the intrinsic anisotropy is 
intimately linked with the Goldreich- Julian limited vacuum inside the pulsar 
magnetosphere where large-scale electric fields can be maintained accelerating 
particles in particular directions. If ambient plasma would shortcut the poten- 
tial drops between open field lines, the emitting particle distributions would 
immediately isotropize. Thus, this elegant way out of the paradox is blocked 
for AGNs which we believe are powered by accretion of high-conductivity 
plasma onto black holes. It is also blocked for GRBs where the energy release 
greatly exceeds the magnetic field energy. 

AgNs: The rapid variability on time scales At < 1 day characterizing 
gamma ray emitting AGNs seems to indicate emission from a scale of the 
order of 10 — lOOrs. A strong rotating magnetosphere cannot be ascribed 
to the black hole itself by virtue of the 'no-hair' theorem, but very likely 
it is carried by plasma accreting onto the black hole. Matter can only spiral 
inward towards the black hole, if its angular momentum is extracted. Angular 
momentum transport could occur through a rotating magnetized wind ejected 
along the rotation axis (Blandford & Payne 1982, Pudritz 1985, Ferreira & 
Pelletier 1993) which is observed in about 10% of all AGN as a powerful 
radio jet. To date, there is no evidence for the competing mechanism of 
radial angular momentum transport in an accretion disk (Kinney 1994). The 
Alfven point on each streamline (the corotation radius) would have to be at 
~ 10 — lOOrs for an ejection/accretion rate ratio of Afg/Afa ~ 1/100. Since 
the released energy is imparted on the small amount of mass ejected through 
the wind, the energy per particle can asymptotically attain very high values, 
i.e. 

r = E/M^c^ ~ L/mIc^ ~ 100?7a ^ 10 (7) 

where L = rjg^Mg^c^ denotes the accretion luminosity with efficiency ?yacc — 0.1 
in the Schwarzschild metric of a non-rotating black hole. In a Kerr metric, the 
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efficiency r]a_ can be signfficantly higher. Relativistic hulk motion is required, 
if a tenuous magnetized wind carries away the angular momentum of the 
accretion flow. Due to the pinching Lorentz force by the self-produced toroidal 
magnetic field, the wind can coUimate to a jet fiow emerging from the central 
engine as a very narrow cone. However, the jets are believed to be very 
different from pulsar beams: The particles and photons in the comoving frame 
of the jet can be assumed to be isotropic due to rapid pitch-angle scattering 
off magnetic field fiuctuations. Beaming of the radiation pattern arises solely 
by bulk relativistic motion with Lorentz factor F ^ f in a direction at an 
angle 6 with the line of sight. From phase space covariance, we obtain the 
observed fiux density (Rybicki & Lightman f979) 

SoU'^) = S'+\S{:y) (8) 

where S{p) denotes the fiux density in the comoving frame and where 6 
denotes the Doppler factor given by 

S = ^ (9) 

T{1- 13 cos 9) 

For the special angle cos 6 = fi we obtain 8 = V . The observed variability 
time scale is given by 

AUs = ^tjS (fO) 

leading to underestimates of the source size r < cAt. Inserting the photon 
density of a source at distance D 

^■y = J^{j)' J S{i^)dlni^ = r(^+«)nobs (11) 
into Eq.(5) yields the optical depth at the source 

T^^^„ = ,5-(^+")T^^,o(obs) (12) 

Typically, a has values from 1/2 to 1 and the optical depth is reduced by 
~ F~^. Thus, for F > Fo ~ 3 the photon-photon pair production optical 
depth decreases below unity. As a matter of fact, the maximum energy of 
blazars detected by EGRET is at least ~ 10 GeV requiring T^-y^o < lO""*" or 
F ~ 8 — 16 in sources like 3C279 which is in accord with the simple estimate 
Eq.(7) for a jet with MJM^ ~ 1/100. 

GrBs: Again referring to Tab.l, the gamma ray absorption problem seems 
to be the worst in the case of gamma rays bursts. Although they emit only 
a short pulse of gamma rays with a typical duration of ~ 1 s, the observed 
fiuence F ~ 10"'' erg cm~^ (Hartmann 1995) implies a luminosity of 

L ~ 10^\D/m Mpc)2(F/10-^ erg cm-^)/{tgrh/l s) erg s"^ (13) 
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Even a halo origin at D ~ 100 kpc reducing the GRB luminosity by six orders 
of magnitude compared to cosmological GRBs does not resolve the paradox. 
For a cosmological distance D ~ 100 Mpc, which is more in line with the 
observed isotropy and apparent brightness distribution of GRBs (Tegmark et 
al. 1995), the total energy release in photons E ~ 10^^ erg is a tantalizing 
1% of the gravitational energy of a collapsing solar mass. In a supernova, the 
energy of the collapsing core is released into a mass shell with M ~ 10 Mq 
provided by the surrounding massive star. Thus, in the latter case the ejecta 
cannot reach relativistic velocities. A more or less naked collapse is required 
to obtain F > 1, e.g. due to the coalescence of two neutron stars (Meszaros & 
Rees 1993a). Then the energy release occurs in the low density region perpen- 
dicular to the orbital plane of the coalescing or accreting neutron star binary 
where only a small baryon pollution is expected. Primarily the energy goes 
into p -\- p pairs and it has been estimated that 1% of the neutrinos annihi- 
late yielding electron/positron pairs, the relativistic 'fireball"*. The Thomson 
thick fireball is highly super-Fddington, i.e. L/Le — 10*"^, and it is therefore 
subject to rapid expansion. During adiabatic expansion, the available energy 
is very soon consumed by the few polluting baryons which store the energy 
as bulk kinetic energy (Meszaros et al. 1993). Radiative losses during these 
stages of the explosion are very small. It is the formation of a coUisionless 
shock when the mass of the swept-up interstellar matter equals the exploding 
mass that non-thermal radiative losses become important. This is analogous 
to supernova remnants, with the difference that the fireball explosion is rel- 
ativistic. The optical depth for a relativistic, spherically expanding shell is 
given by 

T^^^„ = F-(^+")T^^,o(obs) (14) 
and the observed variability time scale by 

AUs = At/F^ (15) 

Hence the criterion T^-y^o < 1 requires a bulk Forentz factor 

To = TT^ = ~ - (16) 

where Afg denotes the conserved baryonic rest mass. The size of the expanded 
fireball is = F^cAt ~ 10*^ cm which is much larger than the initial compact 
object. Its small optical thickness is in accord with the observed non-thermal 
spectra. By contrast, the precursor event is extremely short, optically thick 
and less bright than the optically thin relativistic blast wave (Meszaros & 
Rees 1993b, Meszaros et al. 1994). A GRB should therefore be viewed as a 
relativistic supernova remnant. 

■*For a criticism, see Janka & Ruffert, MPA-Report 909 
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Figure 1: 3C273: Proton-initiated cascade spectrum fitted to the non-thermal 
jet spectrum (Mannheim 1993b). In addition to the jet spectrum, the data 
(Lichti et aL 1995) show thermal emission in the FIR, UV and soft X-ray 
bands presumably emitted by an accretion flow and a surrounding dust torus. 
Above ~100 GeV, the gamma ray spectrum cuts off sharply due to pair 
absorption in the infrared radiation field of the dust torus 

3 Relativistic jets and accreting black holes 
in AGNs 

The most detailed, but still controversial, understanding of the radiation 
mechanisms responsible for gamma ray emission from compact objects has 
been achieved for pulsars (e.g., Cheng et al. 1986, Daugherty & Harding 
1982, Gallant & Arons 1994) and AGNs (e.g., Blandford & Levinson 1995, 
Dermer & Schlickeiser 1993, Henri et al. 1993, Sikora et al. 1994). In this 
Lecture, emphasis is on the author's hadronic emission model for AGNs. The 
non-stellar activity characteristic for AGNs basically comes in two manifes- 
tations, (i) a UV-to-soft-X-ray bump indicating thermal emission and (ii) a 
broad non-thermal continuum extending from radio frequencies to gamma 
rays. The relative strength of the two emission components varies by orders 
of magnitude between members of different classes of AGNs. Radio-quiet 
AGNs show practically only the thermal emission component (plus a thermal 
FIR bump due to reprocessing by heated dust and a weaker X-ray component 
which could be either thermal or non-thermal), whereas radio-loud AGNs of 
the blazar subtype only show the non-thermal component. Fig.l shows the 
mixed thermal/non-thermal spectrum of the quasar 3C273. 
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Radio-LOUD AGNs: Morphologically, at a sufficient angular resolution 
powerful radio emission is related to radio jets. For very large Alfvenic Macli 
numbers Ma ^ 1, the rotating hydromagnetic wind from the accretion flow 
onto a supermassive black hole can emerge from the central parsec of the AGN 
propagating out of the host galaxy and deep into intergalactic space (Camen- 
zind f990). Estimating the kinetic power of radio jets by applying strong 
shock jump conditions to the terminal shocks (hot spots) of jets where they 
hit the intergalactic medium (Roland et al. f988, Meisenheimer et al. f989), 
by computing the radio age (Rawlings & Saunders f99f) or by synchrotron- 
self-Compton limits (Celotti & Fabian f993, Daly et al. f996), generally 
yields values of the order of the Fddington luminosity for black holes in the 
mass range fO® — fO^^Af©. Therefore, it is evident that the jets must emerge 
from the central engine, no other comparable power supply is available. The 
emergence of the jet from the central region has the fascinating implication 
that the plasma from the vicinity of the black hole can be probed on scales 
which are observationally resolvable using very long baseline interferometry 
(VLBI). One of the major results is that knots of enhanced surface bright- 
ness in the radio jets can move at superluminal speeds (e.g., Porcas 1987, 
Begelman et al. 1994). The most natural and robust explanation of this ob- 
servation is that of relativistic motion of the emitting plasma at a small angle 
to the line of sight which can give apparent speeds t^j,app ~ T^^j when projected 
onto the plane of the sky. Another consequence of relativistic bulk motion is 
the beaming of the photon flux into a cone of opening angle l/F, see Fq.(8). 
Consequently, the appearance of a radio source must depend strongly on the 
angle between the radio axis and the line of sight. As a matter of fact, the 
zoo of radio-loud AGNs can be unified to the species of powerful (FRII) radio 
galaxies (seen edge-on) appearing as fiat-spectrum radio quasars (Fig.l) when 
seen pole-on and to the weaker FRI radio galaxies appearing as BL Lacertae 
objects (Fig. 2) for a pole-on orientation (Urry & Padovani 1995). In addition 
to the orientation-dependence of the radio emission, the nuclear emission line 
regions and the UV bump also vary with angle (Antonucci 1993). An obscur- 
ing torus hides the central parsec within ~ 45° from the edge-on orientation. 
In the spectral energy distribution, the importance of the non-thermal con- 
tinuum increases with increasing pole-on orientation. All radio-loud AGNs 
with their jet axes oriented at small angles to the line of sight are summarized 
under the name 'blazars' due to their variable, polarized synchrotron spectra 
from a boosted relativistic jet. The major unanswered questions affecting the 
nature of the gamma ray emission mechanism are: 

(1) What is the composition of the jet plasma (Poynting fiux, pairs, 
solar abundance plasma, enriched plasma, something else)? 

(2) At which distance to the black hole is the emission produced? 
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(3) What is the comoving-frame ratio of magnetic field energy density and 
radiation energy density? 

(4) What is the acceleration mechanism producing the ultrarelativistic par- 
ticles in the comoving frame responsible for the non-thermal radiation? 

In the baryon-induced cascade model for the gamma rays presented in this 
Lecture, it will be assumed that (1) the plasma contains protons and electrons, 
(2) the distance of the non-thermal source in the jet is large enough (~ 
1 lyr) to make the locally produced synchrotron photons the dominant target 
contributing to the gamma ray opacity, (3) the magnetic field energy density 
is assumed to dominate over the radiation energy density and, finally, (4) 
shock acceleration is the acceleration mechanism at high energies. These 
assumptions lead to a consistent physical picture and a predictive model 
(Sect. 4). Let us discuss some arguments in favor of the above assumptions. 

Ad (1) If the jet originates as a hydromagnetic wind ejected by the accre- 
tion fiow, the plasma fiowing along the streamlines is of the same composition 
as the plasma accreted onto the black hole. Protons should be the dominant 
baryonic species in this plasma. Alternatively, if the jet originates as Poynting 
fiux produced by a maximally rotating black hole, the Poynting fiux would 
be converted to a relativistic pair plasma close to the ergosphere by photon- 
photon interactions in the UV/X-ray radiation field of the optically thick 
accretion fiow. As in the case of fireballs, however, the unavoidable polluting 
baryons ultimately receive the released energy during an adiabatic expansion, 
while the pairs would cool and annihilate. Nevertheless, the absence of Fara- 
day depolarization seems to argue in favor of a pure pair plasma in which 
Faraday depolarization is absent (Jones 1988). However, in an ordinary e-p 
jet plasma which is relativistically hot, with temperatures ~ 100 MeV, Fara- 
day depolarization would also be absent. Such temperatures are expected 
behind shock waves propagating down the tenuous jet fiow (Blandford & 
McKee 1976). Runaway pair production by coUisional processes does not oc- 
cur, since the density in the jet plasma is too low (particle- wave interactions 
dominate). 

Ad (2) Dissipation is certainly important at the base of the jet where an 
entropy generating demon loads particles onto the field lines. The heated jet 
base would act as a comptonizing screen through which the UV photons from 
the optically thick accretion fiow have to pass. Unsaturated comptonization 
of the UV photons by the hot jet base produces power law soft X-ray excesses 
which are commonly observed in AGNs (Walter & Fink 1993). For the par- 
ticle energy distributions to develop power law tails reaching ultrarelativistic 
energies, it is necessary that the density in the jet becomes very low and that 
collisionless shocks form. While the formation of shocks due to magnetic fiare 
activity at the base of the jet seems to be unescapable, the rather high density 
n ~ 10^° — 10^^ cm~"^ and consequently large Coulomb coUisional rate tends 
to thermalize the injected energy (tcoui < n/*^ ^ ^ 1^^° cm~"^ at T = 10^ K 
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Figure 2: Mrk421: Proton blazar fits to tlie liiglily variable broadband con- 
tinuum obtained by changing the angle between the velocity of the emitting 
plasma in the jet and the line of sight by a few degrees (shocks traveling in dif- 
ferent directions?). Note that Mrk421 is detected over 18 orders of magnitude 
in frequency 

and Tj = lOOrs, Mannheim & Schlickeiser 1994). In fact, flare-induced shock 
dissipation could be the major heating mechanism of the jet base presum- 
ably responsible for soft X-ray emission. Therefore, the most likely site of 
the gamma ray emission region in the jet is far away from the base of the 
jet when n <^ 10^° cm~"^. Shocks would be generated behind the magnetic 
nozzle provided by the self-confining toroidal field (van Putten 1995) when 
the jet expands into the steep pressure gradient of the interstellar medium of 
the host galaxy. The distance of this region to the central black hole could 
be very large, i.e. r ^ 1 la. 

Ad (3) The most direct measurement of the ratio a = u^/ub comes from 
the hot spots of Cyg A (Harris et al. 1994) yielding a < 10~^, cf. Fig. 3. 
The idea behind the method is that X-ray observations limit the comoving- 
frame ratio of the photon-to-magnetic energy density, since the observed syn- 
chrotron emission from the jets implies synchrotron-self-Compton X-rays due 
to the Compton scattering of synchrotron photons with energy €syn off the 
same relativistic electrons with Lorentz factor 7e which produced them in the 
first place. In the Thomson scattering regime (esyn7e < mgC^), the character- 
istic photon energy of inverse- Compton scattered photons is given by 

eic ~ esyn7e (17) 
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Figure 3: Cyg A: Proton-initiated cascade fits to the hot spot spectra. The 
radio-to-optical spectrum is synchrotron emission by accelerated electrons. 
The X-ray-to-gamma-ray spectrum is unsaturated synchrotron cascade emis- 
sion by accelerated protons assuming Up = 40Me and i? = 3 X lO""* G 



which easily reaches the X-ray range. The SSC X-ray luminosity can be 
approximated by 

-^ssc ~ —Lsy-a (18) 
Ub 

In a stationary state, the photon energy density should not exceed the mag- 
netic field energy density, since the particles are required to scatter ofi' the 
magnetic field fiuctuations usb ~ ub to be accelerated by a statistical Fermi 
(e.g., Drury 1983) or drift mechanism (Webb et al. 1983, Begelman & Kirk; 
1990). Thus, if u-^ > ub the emitting particles would escape from the mag- 
netic confinement. Moreover, the back-reaction of the inverse- Compton cool- 
ing on the electron energy distribution would act to equilibrate the energy 
losses over a broad energy band leading to continua smoothly joining the 
optical and X-ray ranges by an a = 1 power law which is apparently not 
the case for the Cyg A and other hot spots. Assuming now u-^ oc and a 
tangled magnetic field B oc (adiabatic expansion in a conical jet implies 
i?_L oc and B\\ oc r~^) leads to a = constans implying that synchrotron 
emission is the dominant cooling mechanism for the relativistic electrons in 
the jet. This is in accord with the fact that the beamed hard X-ray lumi- 
nosity from nuclear jets typically does not exceed the synchrotron luminosity 
Fq.(18). Nevertheless, it could be possible that the external radiation field 
energy density exceeds the magnetic field energy density close enough to the 
thermal source of UV photons. This would not be the case in angular mo- 
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mentum carrying jets where an initial equipartition = ub would lead to 
an ever-increasing domination of the magnetic field in the jet with increas- 
ing distance^. Thus models in which the dominant energy loss is not due to 
synchrotron emission imply the assumption that the jet magnetic field is ini- 
tially below equipartition and that it is increased by dynamo action in the hot 
spots. The importance of the rotating jet as an angular momentum transport 
structure is then diminished accordingly. 

Ad (4) Shock acceleration theory has been applied successfully to hot 
spot spectra at the kpc scale and the local spectra produced by knots in the 
jet at smaller scales (VLA, VLBA, VLBI) seem to be quite similar. Moreover, 
there is a detailed correspondence of shock-in-jet models with the observed 
variability (Hughes et al. 1989, Quian 1991) and from numerical simulations 
of supersonic jets crossing steep pressure gradients shock formation seems 
unavoidable. Periodic shock structures can be formed by magnetic nozzles 
(van Putten 1996). Moreover, shock acceleration is known to work from in 
situ experiments in the solar system and it is the leading contender explaining 
the acceleration of cosmic rays at supernova remnants. The major drawback 
is the lack of a detailed understanding of how relativistic and oblique shocks 
accelerate particles. A large computational and theoretical effort is necessary 
to improve the situation. 

Radio-QUIET AGNs: An explanation of the absence of strong gamma 
ray emission from radio-quiet AGNs (Maisack et al. 1995) is a substantial 
requirement for any model explaining gamma rays from radio-loud AGNs. In 
radio-quiet quasars and Seyfert galaxies, the thermal big blue bump domi- 
nates the bolometric luminosity by far (Sanders et al. 1989). Powerful radio 
jets are generally not present in these objects, but faint bipolar radio mor- 
phologies seem to be quite frequent. As in the case of radio-louds, accretion 
onto a black hole and angular momentum transport by a rotating MHD wind 
constitutes a viable explanation. The competing idea of a thin accretion disk 
where angular momentum is transported radially outwards by viscous friction 
between turbulent eddies is in confiict with observations (Kinney 1994). E.g., 
the optical emission should come from much larger radii than the UV and soft 
X-ray emission. Nevertheless, fiux variations occur quasi-simultaneously from 
the optical to X-ray range (Clavel et al. 1992). Coupled accretion/ejection 
structures could ameliorate the problem in many ways (Mannheim 1995b). 
On the other hand, the absence of radio (and gamma ray) emission from radio- 
quiets is often invoked as counter-evidence against a nuclear wind. However, 
if the Alfvenic Mach number of the rotating magnetized wind is ~ 1, Kelvin- 
Helmholtz instabilities grow rapidly over a length scale of order ~ lOrs. As 
a consequence, the vertical wind with large opening angle ~ I/Ma becomes 

®For comparison, a thermal radiation field of the same power as the jet in a typical 
FR II radio galaxy with L ~ 10'*® erg s~^ would have an energy density of L/iAirD'^c) ~ 
3 X 10~^^ erg cm~^ at a distance of _D ~ 100 kpc, whereas the magnetic field B ~ 10"'* G 
has B^/lSir) ~ 4 x 10~^° erg cm~^ in hot spots. 
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turbulent (Mannheim 1995b). In this way, the formation of strong shocks 
propagating down the jet producing powerful radio emission is avoided ex- 
plaining 'radio-quietness'. The wind in its outer parts is colder and only par- 
tially ionized providing a warm absorber seen in the X-ray spectra of many 
AGNs. As the formation of the hydromagnetic wind is certainly accompanied 
by dissipative processes, the base of the tenuous wind close to the putative 
black hole must be heated to very high temperatures, e.g. by magnetic re- 
connection. Any line of sight to the UV-emitting optically thick plasma of 
the accretion flow must pass through this hot wind base, thus producing soft 
X-ray power law spectra by unsaturated Comptonization of the UV photons 
- there is no fundamental difference between radio-quiets and radio-louds in 
this respect (the soft X-ray spectra in radio-louds should be slightly flatter). 
Rapid correlated variations on a dynamical time scale occur naturally in this 
scenario (e.g., due to coronal mass ejections). Following the arguments in 
Sect. 2, radio-quiet AGNs should represent disk/jet systems where the ejected 
mass loss rate larger than in sources with relativistic jets and consequently 
their jets are non-relativistic. Any non-thermal emission associated with par- 
ticle acceleration in the wind (e.g., stochastic acceleration Henri & Pelletier 
1991, Dermer et al. 1995, Mannheim 1995b) occurs well within the gamma 
ray photosphere of the thermal UV/X-ray radiation field (t^^ > Ij; since 



assuming = 1, Zx = 0.1e_iivE and the scattered isotropic radiation field 
-^x,iso = 0.0lT_2ivx. Whereas hard X-rays are refiected by Thomson scat- 
tering off the cold accretion fiow, the gamma rays are absorbed heating the 
fiow. Therefore, radio quiet AGNs cannot be powerful gamma ray emitters^. 
The absence of a diffuse gamma ray bump with a power comparable to the 
30 keV X-ray bump in the diffuse background radiation field supports these 
conclusions. 

4 Baryon-initiated cascades 

The spectra of radio jets are characterized by a polarized radio-to-X-ray fea- 
tureless continuum which is commonly interpreted as synchrotron emission 
by relativistic electrons with energies of maximally ~ 100 GeV. Typically, 
the spectral slopes of local emission regions in the jet (knots, hot spots) are 
a ~ 0.5 (5'j, oc in the radio-to-submm range and a ~ 1 in the infrared. 

Further steepening is observed in the the optical-to-soft-X-ray range. At 

^There is a possibility that soft gamma rays of energy e ~ few MeV are emitted with 
an anisotropic distribution preferentially towards the disk. In this case, we expect weak 
gamma ray emission from Seyfert 2 galaxies which we believe are AGNs seen near the plane 
of the accretion flow. 




(19) 



14 



Unsaturated Proton-Initiated Cascades 
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Figure 4: Sketch of the PIC mechanism: Shock-accelerated electrons produce 
radio-to-X-ray synchrotron radiation (target). Shock-accelerated protons cool 
in the synchrotron target photons and magnetic field producing pions, pairs 
and synchrotron emission (for the sake of demonstration, all three secondary 
luminosities have been made equal). The neutral pions decay producing ul- 
trahigh energy gamma rays [top panel). Electromagnetic power injected into 
the acceleration region above the energy E* where t^^ = 1 is absorbed cre- 
ating relativistic pairs. These pairs initiate the first generation of cascade 
gamma rays via synchrotron emission [middle panel) and the cycle repeats. 
In the above example, already the second generation of gamma rays emerges 
almost completely below E* terminating the cascade [bottom panel). The re- 
sulting total spectrum is given by the superposition of the individual absorbed 
cascade generations and the primary electron synchrotron spectrum. 
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hard X-rays, the spectra seem to rise again possibly indicating synchrotron- 
sell-Compton emission. The origin ol hard X-rays and gamma rays cannot 
be resolved spatially, however, their association with jets seems to be beyond 
doubt. In greatest detail, the radio-to-X-ray spectra have been studied for 
hot spots in the jets of FR II radio galaxies at a scale of ~kpc. By two inde- 
pendent methods the magnetic field strengths in hot spots could be inferred 
yielding i?hs ~ 3 X lO""* G (Meisenheimer et al. 1989). Closer towards the 
core, the spectra of knots in the jet show a low-frequency synchrotron-self- 
absorption turnover increasing in frequency with decreasing distance from 
the core. Apparently, regions of enhanced surface brightness density along 
the jet represent scaled-down versions of hot spot spectra in regions of larger 
magnetic field (Bj ~ 0.1 — 1 G). Knots at very small scales ~ 0.1 — 1 pc 
superimpose to a more or less conical emission region. As a result of the 
superposition of increasingly self-absorbed spectra, the nuclear jet obtains a 
fiat radio spectrum (Kellermann & Pauliny-Toth 1981). VLBI observations of 
superluminal motion in nuclear jets indicate relativistic bulk motion, whereas 
the more distant parts of the jet (hot spots) are very likely at least moderately 
relativistic resulting from a general slowing down of the jet. Hence it follows 
that in the relativistic part of the jet the hot spot spectra are boosted by the 
Doppler factor (5 ~ 3 — 10 in frequency and by 6"^ in luminosity. Therefore, the 
nuclear jets appear as extremely bright and variable blazars when viewed at 
a small angle of the jet axis to the line of sight. On the other hand, when the 
observer is outside of the boosting cone with opening angle l/F, the nuclear 
jet is greatly diminished in brightness and the otherwise washed out more 
distant and isotropically emitting parts of the jet (hot spots, lobes) become 
visible. In spite of their brightness, blazars are very difficult to diagnose be- 
cause of the relativistic projection effects, the lack of spatial resolution at 
a high dynamical range and the complex environment in the center of the 
active galaxy. Therefore, the best laboratory to study the non-thermal radi- 
ation processes in jets are the isolated hot spots of FR II radio galaxies, such 
as Cyg A. 

What is the origin of the relativistic synchrotron-emitting electrons? For 
radio jets which reach out to distances of ~Mpc, it is impossible that the 
relativistic electrons (or positrons) originate in the active nucleus. Before 
they arrive at the hot spots, they would loose their energy by inverse- Compton 
scattering off the 3K-microwave background. Therefore, in situ acceleration 
in the hot spots is required. The observed power law energy distributions 
dn^ oc ^f~'^d^fQ argue for a statistical process. Diffusive (e.g., Drury 1983) 
or drift transport (Begelman & Kirk 1990) of relativistic particles across a 
strong coUisionless shock are the leading contenders among the possible in 
situ acceleration mechanisms. A major drawback of the theory is the lack of 
an efficient pre-acceleration mechanism. Resonant scattering off Alfven waves 
requires that the particle momenta exceed a threshold momentum which is 
already relativistic in the case of electrons. Pre-acceleration of electrons out 
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Figure 5: Sketch of the blazar spectrum due to emission from varius nuclear 
jet radii. Injection of a shock at small radii propagating along the jet increases 
the acceleration efficiency and maximum energies producing a 'flare'. The 
parameter ^ denotes the ratio of the proton maximum energy to the limit in 
Eq.(26). For fresh shocks, ~ 1 while nonlinear shock evolution decreases 
the acceleration efficiency yielding < 1 



of the thermal pool could be achieved by resonant interaction with Langmuir 
and Whistler waves, or by turbulent reconnection of magnetic field lines. For 
protons, the injection problem is less problematic, since a rather large number 
of protons from the Maxwell tail of the thermal distribution is available above 
the Alfvenic threshold (Malkov & Volk 1995). Moreover, the Larmor radius 
of a gyrating proton is much larger than that of an electron increasing the 
acceleration efficiency, since the shock front is likely to be broadened to a 
nonlinear structure. The Larmor radius sets a typical scale for pitch-angle 
scattering, so traversing the shock on a single scattering (which is essential 
to gain energy from the Lorentz transformation between the upstream and 
downstream side of the shock) should be easier for protons than for electrons. 
In the local cosmic ray gas, the energy density of relativistic protons above 
GeV exceeds that of the relativistic electrons by a factor of ~ 100, although 
this does not prove that proton acceleration is generally more efficient at high 
energies. A first-principle prediction of the p/e energy density ratio rj = Up/ue 
is not possible at the current status of the theory. 

Historically, protons have generally been neglected as a radiative agent. 
Their main energy loss at high energies is pion production by collisions with 
matter as indicated by the 7r° bump at 70 MeV in the gamma ray emission 
from the galactic disk due to the interaction of cosmic rays with the interstel- 
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lar matter. However, in radio jets the matter density is extremely low, viz. 
~ 10"'* cm""^ at r ~ 1 kpc. The p-p pion luminosity for a Lsyn ~ 10"*^ erg 
jet amounts only to the negligible value 

Xpp ~ ^ ~ vLsyJ-f' ~ lO-'vLsy. (20) 
tpp tpp 

The negligible energy losses of relativistic protons imply that they must reach 
very high energies in a statistical process balancing escape and acceleration. 
In fact, the energies reach up to values where the acceleration time scale 
equals the energy loss time scale (unless the Larmor radius would exceed the 
jet radius). As pointed out by Biermann & Strittmatter (1987) and Sikora et 
al. (1987), protons should therefore reach energies where the acceleration is 
halted by photo-production of pions and pairs. The processes 

p -\- ^ —i- Tr° -\- p , Tr~^ -\- n (21) 

p + -f^e+ + e-+p (22) 

have effective inelastic cross sections cr^ ~ m^/mpCtfCTpp ~ 10~^® cm^ and 
(Te± ~ 2me/mpaf(TT ~ 10""^° cm^ (Blumenthal 1970). The dimensionless 
threshold energies for a head-on collision with a low energy synchrotron pho- 
ton of energy e are 

m^c2(l + f^) 

^ 2mp / 

7p,th,7r = — (2dj 

and 

2meC^ 

7p,th,e± = (24) 

Typical values for optical target photons e ~ 1 eV are 7p,th,7r ~ 10^ and 
7p,th,e± ~ 10^. In a power law target radiation field with a = 1 (IR-to- 
optical), pair and pion production contribute equally to the proton cooling. 
The proton maximum energy is given by equal acceleration and energy loss 
time scales. Since the shortest acceleration time scale is obtained for Bohm 
diffusion (tacc ^ '"l/c), we arrive at the estimate 

^'■""S2xl0"(j^^^)"%a (25) 

Additionally, we require that the Larmor radius does not exceed the jet radius 
(Bell 1978) obtaining 

Therefore, proton acceleration in hot spots (a ~ 0.01) is bounded to values of 
order 7max ~ 3 X 10^^ by geometrical constraints. In the comoving frame of 
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nuclear jets (i? ~ 1 G, r ~ 0.1 pc) we obtain 7inax ~ 3 X 10^°(a/0.01)~^/^ for 
Bohm diffusion and the Larmor radii of the most energetic particles remain 
an order of magnitude below the jet radius. 

The proton secondary pairs and pions are the origin of non-thermal high 
energy radiation from jets. The pairs cool instantaneously via synchrotron 
radiation (inverse- Compton radiation is down by a factor a ~ 0.1 — 0.01), the 
pions decay into gamma rays via 7r° — )• 7 -|- 7 and produce further positrons 
via TT e"*" -|- 3p. We have seen in Sect. 2 that the optical depth for gamma 
rays of energy E is proportional to E". Since E reaches very large values for 
gamma rays of secondary origin, the optical depth inevitably exceeds unity. 
The gamma rays are absorbed within the radio jet producing pairs, which 
themselves produce further gamma rays etc. The cycles terminate when the 
gamma ray energy of the n-th generation decreases below the energy E* where 
T^^[E*) = 1. The process is sketched in Fig. 4 (Mannheim et al. 1991). Note 
that the thermal pair yield, and hence the flux of the annihilation line, is 
very low, since unsaturated cascades produce only highly relativistic pairs 
(Svensson 1987). 

In the most simple model of the cascade, local power law injection of sec- 
ondary pairs and gamma rays as well as primary electrons into a spherical 
region is assumed. The emerging photon spectrum corresponds to a station- 
ary solution of the coupled set of kinetic equations for the pair energy dis- 
tribution n(7,t) and the photon energy distribution m[x^t) at dimensionless 
energy x = e/nicC^ given by 



where the qi denote sources and sinks due to pair production and injection and 
the other terms denote energy losses (emissivities, resp.) due to synchrotron 
and Compton radiation. Diffusive escape from the acceleration region is ac- 
counted for by an escape time which is related to the lowest eigenvalue of the 
spatial diffusion equation. Acceleration would be encompassed by including 
diffusion terms and by solving the kinetic equations in the shock geometry. In 
Eqs.(27) and (28), the acceleration process is hidden in the power law injec- 
tion terms. However, in the high-energy limit (X-rays, gamma rays) diffusion 
is not the dominant process for the secondary pairs. 




(27) 




7, esc 
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Figure 6: The cosmic diffuse neutrino baclcground produced by AGNs. The 
diagonai fine shows the atmospheric baclcground due to iarge-angie //'s from 
cosmic ray interactions (Gaisser et ai. f 993). Experimentai iimits by the Fre- 
jus and FAS-TOP coiiaborations are shown for comparison (see Mannheim 
f995a for references). The PeV limit comes from the non-detection of con- 
tained events due to the Glashow- Weinberg resonance assuming an equal 
mixture of electron and muon neutrinos 



In Fig. 3 the model is applied to the hot spots of Cyg A where X-ray 
emission was detected with the High Resolution Imager of ROSAT (Harris et 
al. 1994). Synchrotron-self-Compton emission is neglected using a magnetic 
field value of 3 X 10"'* G and r/ = 40. Neglecting proton-initiated-cascade 
(PIC) emission would yield a magnetic field of 10""* G for i] = 0. The spectral 
slope of the SSC emission would be close to a = 1, whereas the PIC emission 
has a = 0.7. The sensitivity of current gamma ray experiments is insufficient 
to test the model predictions at very high energies. 

To model the entire blazar spectrum with the simple spherical proton- 
initiated-cascade model, we proceed as follows. The total primary electron 
synchrotron spectrum is modeled as a superposition of a large number of in- 
dividual knots in a conical jet with the magnetic field scaling as i? oc 
(Blandford & Konigl 1979). For proton-initiated cascades, only target pho- 
tons in the infrared-to-optical regime are important due to the threshold en- 
ergy Fqs.(23),(24). The largest radiation compactness of infrared-to-optical 
photons comes from the jet radius rb where the spectral steepening of the local 
synchrotron spectrum due to cooling across the expanding emission volume 
occurs at the same frequency as the synchrotron-self-absorption (Blandford 
& Konigl 1979). Cascade emission from larger radii emerges at higher en- 
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ergies (lower value of t^^) and has a lower flux (less target radiation), see 
Fig. 5. At r < Tb the synchrotron efficiency of the jet is greatly reduced 
due to self-absorption. Therefore, modeling cascade emission by a Doppler- 
boosted spherical blob at r = rb is a fair approximation. A flare due to the 
propagation of a strong shock along the conical jet would reach its maximum 
luminosity at this stage. All relevant physical parameters of the jet except for 
the proton/electron ratio rj can be obtained from their radio-to-X-ray spec- 
tra and morphology. Together with gamma ray measurements, the 'proton 
blazar' model allows for a complete diagnosis of the jet plasma. 

Results from fitting the proton blazar model to a number of nearby blazar 
spectra are consistent with a more or less constant proton-to-electron energy 
density ratio rj ~ 100 (Mannheim et al. 1995b). The differences in gamma 
ray luminosities seem to be largely due to differences in the proton maximum 
energies. Variability can occur on very short time scales for fiattened Mach 
disk like geometries, inhomogeneities in the jet fiow (Quian 1991), or opti- 
cal depth effects of the gamma ray photosphere. Gamma rays and optical 
emission should vary in a correlated way. Short time lags (photo-production 
cooling time scale, proton acceleration time scale) are possible. The gamma 
ray variations have a larger amplitude than the target photon variations due 
to the quadratic effect that the gamma ray fiux is dependent on 7p,max and on 
7e,max, whcrcas the radio-to-X-ray synchrotron emission is only proportional 
to 7e,max- Mouitoring of 3C279 is consistent with these predictions (Maraschi 
et al. 1994). Due to the often sparse frequency coverage of the measurements, 
there is still considerable ambiguity in the diagnostic method. In addition, 
the blazar spectra are highly variable demanding contemporaneous data sets. 
It is therefore a prime goal to coordinate contemporaneous observational cam- 
paigns on blazars in the radio, mm, submm, FIR, NIR, optical, UV, X-ray 
and gamma ray ranges. 

5 High-energy neutrinos and their detection 

Decay of the charged pions produces high energy neutrinos 

TT^ ^ fi^ + Vf, ^ + Ve + Vf, + Vf, (29) 

Typically, the four light particles in the final state each receive an energy 
of Ep/20. Considering the maximum energies along the jet Fq.(25), we can 
expect a neutrino fiux peaking in the range 10® — 10^° GeV. The neutrino to 
gamma-ray fiux ratio is determined solely by decay kinematics and the relative 
production cross sections of charged and neutral pions (given by Clebsch- 
Gordan coefficients). As a result, neutrino fiux predictions can be scaled to 
the known photon fiuxes of the sources. In contrast to the photons, however, 
which suffer cascading, the neutrinos are always optically thin. Therefore 
the neutrino spectrum is the original production spectrum. Fig. 6 shows the 
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Table 2: Muon events per 10^ per year and per steradian at various 
threshold energies for the diffuse AGN backgrounds shown in Fig. 6 

Ej, [TeV] radio-loud radio-quiet 

f l5 887 

fO 3.0 836 

fO^ 0.6 458 

fO^ O.f 



predicted neutrino spectra assuming that the diffuse X-ray background is 
made by radio-quiet AGNs and the diffuse gamma ray background by radio- 
loud AGNs. In the case of radio-quiets, it is assumed that the X-rays are non- 
thermal and associated to (anisotropic) pair cascades in low Mach number 
jets. Stochastically accelerated protons in the turbulent jets reach maximum 
energies of ~ 10'' GeV limited by threshold pion production Eq.(23) in the 
UV field of the accretion fiow (Colgate 1983). 

The spectra thus obtained for individual AGNs are integrated over a dis- 
tribution of AGNs in redshift. The limit from the Frejus proton decay exper- 
iment seems to rule out that proton-initiated cascades produce the X-rays in 
radio-quiets (on the assumption that the diffuse X-rays are due to AGNs). 
This is in accord with the fact that the general indication of non-thermal pro- 
cesses in radio-quiet AGNs is very poor. On the other hand, the prediction 
of neutrinos from radio-loud AGNs is a straightforward consequence of shock 
acceleration. Current fiux limits are still off the predictions by orders of mag- 
nitude. However, the linearly rising neutrino-matter interaction cross section 
(up to the energy where the W-boson mass becomes important in the scatter- 
ing propagator) and the long mean free path of the muons make it possible to 
build underwater or underice detectors providing sufficient sensitivity to dis- 
cover the neutrinos. The H2O acts as a Cerenkov medium for the relativistic 
muons resulting from charged current interactions between the extraterres- 
trial neutrinos with matter in the earth. The Cerenkov light is emitted in a 
narrow cone along the direction of the initial neutrino. It can be detected and 
the track can be reconstructed with conventional photomultiplier techniques 
(Gaisser et al. 1994). A class of such neutrino telescopes aiming at 10^ m^ 
effective area are under construction (see references in Stenger et al. 1992). 
In Tab. 2 the expected annual event rates per steradian arising from AGNs 
are given. Scaling to the corresponding diffuse photon fiuxes leaves relatively 
little room for speculation regarding the energy fiux in the diffuse neutrino 
radiation if the photon fiuxes indeed result from proton-initiated cascades. 
The putative AGN fiuxes provide the most important background for direct 
searches of weak-interaction decay channels of dark matter candidates. 
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Figure 7: The gamma ray horizon due to photon-photon pair absorption for 
two values of the Hubble constant [upper curve: Ho = 100 km s~^ Mpc~^, 
lower curve: Ho = 50 km s~^ Mpc~^ 

6 The gamma ray horizon 

High energy gamma ray photons from extragalactic sources must propagate 
through a diffuse soup of low energy photons produced by other galaxies and 
by the big bang before they reach the observer (Fig. 8). If the gamma ray 
photon energy is high enough, i.e. larger than E = l/(e(l -\- /0.5 eV) TeV, 
there is a probability of 1 ] that the photon will be absorbed by 

pair creation 7 + 7 — )• e"*" + e~ (cf. Sect. 2). Since the pair creation optical 
depth rises with energy as t^^ oc E" for a power law distribution of target 
photons, cosmic absorption causes a steep exponential turnover above the 
energy E* where t^^[E*) = 1 due to external absorption (Fig. 9). Since t^^ ~ 
ny^cr^^D^ a measurement of E* entails information on the product ny^D oc 
n-h/ Ho- If the diffuse infrared background would be known from observations, 
one could therefore infer the Hubble constant at an unprecedented distance 
scale where the Hubble flow is faster than any possible mass flows due to 
inhomogeneous mass distributions (De Jager et al. 1994). On the other hand, 
since a direct measurement of ny^ is difficult due to zodiacal and cirrus light 
pollution, one can infer the photon background instead assuming a given value 
of the Hubble constant (Stecker et al. 1992). The photon background density 
nb is a strong function of galaxy evolution: Farly galaxy formation (cold dark 
matter) implies very large values of nb, whereas late galaxy formation (cold + 
hot dark matter) implies low values of ny^. For the computation of the optical 
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Figure 8: The cosmic background radiation from microwaves (3K radiation) 
over FIR photons emitted by dust to UV photons emitted by early galaxies 



depth we use the geodesic radial displacement function 

dl 



dz Ho 



[{1 + z)E{z)]-' 



(30) 



where E[z) is given by Fq.(13.3) in Peebles (1993). For a cosmological model 
with = 1 and A = the function E[z) simplifies to (1 + zY^'^. Hence we 
obtain the optical depth 



c 



dz— / dji 

dz J-i 2 



oo 



Jenb(e)(l + 2;) cr^^{E^€^ fi^ z) 

2 rp /'CO 



/ dz{l -\- zY'"^ / dx— / den\y{e)a^^{E^e^x — l^z) (31) 
Jo Jo 2 JtiY, 



for a non-evolving present-day background density nb, i.e. n'^{z, e')de' = (1 + 
zYnh{^)de where the dash indicates comoving-frame quantities. The shape of 
the present-day diffuse background density nb(e) is obtained by averaging over 
various galaxy formation models presented by MacMinn & Primack (1995). 
We multiplied this shape by a small factor to obtain agreement with the 
background density (including the contribution from the 2.7 K microwave 
background) estimated by Beichman & Helou (1992) in the FIR (for a modest 
galaxy luminosity or density evolution 7 = 2), and by Madau & Phinney 
(1995) in the NIR through UV range (Fig. 8). The predicted background 
densities depend sensitively on which density or luminosity evolution and 
maximum redshift are assumed and can vary by an order of magnitude. Hence 
it is very important to tighten existing constraints on the actual background 
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Figure 9: Templates of absorbed gamma ray spectra assuming for simplicity 
pF^ = const, and no intrinsic cutoff [Ho = 50 km s~^ Mpc~^) 



density (Biller et al. f995) by obtaining better galaxy luminosity functions 
from deep galaxy surveys or by direct measurements of blazar gamma ray 
(external!) turnover detections. 

We numerically integrate the optical depth function and solve for the 
gamma ray horizon t^^[E^z) = 1. Results are shown in Fig. 7 for two values 
of the Hubble constant. Obviously, to measure the diffuse photon background 
from the era of galaxy formation requires a measurement of gamma ray ab- 
sorption over a wide range of gamma ray energies from fO GeV (z ~ 2) to 
30 TeV (z ~ O.Of). The large number of gamma ray emitting AGNs discov- 
ered by FGRFT provide the beacons at the gamma ray horizon which can 
be used for this challenging program. The minimum and maximum redshifts 
of these blazars are z^i^ = 0.0308 (Mrk42f , seen at TeV by two independent 
groups) and z^a_x = 2.f72 (S5 0836-|-7f0). For the important gamma ray 
range from fO-fOO GeV, the silicon strip detector GLAST would be of great 
importance. In the energy range from O.I-IO TeV ground-based Cerenkov 
telescopes and in the 10-100 TeV range scintillator arrays such as HFGRA 
should be dedicated to measure the gamma ray horizon Fig. 7 using blazars 
at different redshifts. 

Another application of gamma ray absorption measurements is to con- 
strain the distance of GRBs (Alexandreas et al. 1994, Hurley 1996). If the 
bursts are distributed to redshifts of z^a_x ~ 2 as indicated by their apparent 
brightness distribution (Cohen & Piran 1994), and if they have high energy 
power law tails as suggested by the detection of GRBs in the field of view of 
FGRFT (Hurley et al. 1994, Dingus 1995), the number of GRBs at a given 
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Figure 10: The number of gamma ray bursts at high photon energies. An 
excess of bursts over the predicted distribution would indicate a local origin 
in a galactic halo, a lack of bursts would indicate that intrinsic cutoffs are 
important. Upper curve: Ho = 100 km s~^ Mpc~^. Lower curve: Ho = 
50 km s~^ Mpc~^. 

gamma ray threshold energy Eo is greatly reduced below the number visible 
at BATSE energies (~ 100 keV). All GRBs at distances greater than the 
gamma ray horizon are pair-absorbed by the intervening low energy photon 
soup. Thus, if the curve shown in Fig. 10 can be measured by gamma astro- 
nomical methods, this would render a local origin of bursts in an extended 
halo highly unlikely. If the number of bursts decreases more rapidly than 
shown in Fig. 10, this would argue for cutoffs intrinsic to the burst spectra 
and would not allow for a conclusion regarding the distance distribution of 
bursts. The third possibility is that the observed number of bursts at high 
energies is larger than shown in Fig. 10. This would unequivocally rule out 
cosmological models of GRBs. F.g., a single detection of a burst at scintillator 
energies 50 — 100 TeV would be sufficient to rule out a cosmological origin. 

7 Conclusions 

Relativistic winds and the formation of coUisionless shocks play a key role in 
the production of the extremely bright non-thermal gamma ray emission from 
many compact astrophysical objects. It may be conjectured that relativistic 
coUisionless shocks are highly efficient in producing ultrarelativistic particles. 
A major theoretical effort to model particle acceleration at relativistic shocks 
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is mandatory. The specific radiation mechanism responsible for the gamma 
rays is still controversial and depends on which particle species is accelerated 
to which maximum energy. One possibility is that protons reach ultrahigh- 
energies due to their negligible energy losses at low energies. The main proton 
energy loss in coUisionless shocks is the photo-production of secondaries in 
the high density synchrotron radiation field produced by the accelerated elec- 
trons. Unsaturated synchrotron cascades induced by the protons at energies 
10® — 10^^ GeV are characterized by a generic spectrum rising in the X-ray 
range and reaching energies of up to ~ 10 TeV. The spectra match the ob- 
served broadband continua of blazars and may indicate shocks in relativistic 
jets. Contemporaneous multi-frequency observations of gamma ray emitting 
blazars from radio frequencies to TeV photon energies will improve our under- 
standing of particle acceleration and will probe the composition of the plasma 
emerging from the vicinity of the putative black holes residing in AGNs. 

The highest proton energies of order ~ 3 X 10^^ GeV are most likely reached 
in hot spots of FR II radio galaxies implying (i) a non-thermal X-ray spectrum 
fiatter and stronger than synchrotron-self-Compton radiation (Mannheim et 
al. 1991) and (ii) a fiux of escaping protons which could explain the local 
ultrahigh-energy cosmic rays (Rachen & Biermann 1993). The recently dis- 
covered non-thermal X-rays from the hot spots of Cyg A and the indication of 
an excess of ultrahigh-energy cosmic rays from the supergalactic plane (Stanev 
et al. 1995) are in agreement with the predictions. However, a synchrotron- 
self-Compton origin of the X-rays is ruled out only if i?hs > 2 X 10""* G and 
the statistical significance of the ultrahigh-energy cosmic ray/supergalactic 
plane association is not yet fully compelling. An unambiguous signature of 
proton acceleration in AGNs is the emission of ultrahigh-energy neutrinos. 
Neutrino detection with underwater and underice detectors currently under 
construction is feasible. The predicted astrophysical sources of neutrinos pro- 
vide the most important background for direct searches of supersymmetric 
dark matter candidates. 

Irrespective of the nature of the gamma rays from AGNs and GRBs, their 
absorption by intervening diffuse radiation fields can be used as a plain probe 
of the cosmological distance scale or the density of the diffuse far-infrared- 
to-ultraviolet photons. A program to constrain models of galaxy formation 
requires blazar observations in the energy range 10 GeV to 10 TeV with 
satellite gamma ray detectors such as GLAST and ground-based detectors 
such as HEGRA and Whipple. 



References 

Alexandreas, D.E., Allen, G.E., Berley, D., et al., 1994, ApJ 426, LI 
Antonucci, R., 1993, ARA&A 31, 473 

Axford, W.J., 1981, Proc. Int. Cosmic Ray Conf. (Paris) 12, 155 



27 



Begelman, M.C., Kirk, J.G., 1990, ApJ 353, 66 

Begelman, M.C., Rees, M.J., Sikora, M., 1994, ApJ 429, L57 

Beichman, C.A., Helou, G., 1991, ApJ 370, LI 

Bell, A.R., 1978, MNRAS 182, 443 

Biermann, L., 1951, Zeitschrift f. Astrophysik 29, 274 

Biermann, P.L., Strittmatter, P.A., 1987, ApJ. 322, 643 

BiUer, S.D., Akerlof, C.W., Buckley, J., et al., 1995, ApJ 445, 227 

Blandford, R.D., McKee, C.F., 1976, Phys. Fluids 19, 1130 

Blandford, R.D., K6nigl, A., 1979, ApJ 232, 34 

Blandford, R.D., Levinson, A., 1995, ApJ 441, 79 

Blandford, R.D., Payne, D.G., 1982, MNRAS 199, 883 

Blumenthal, G.R., 1970, Phys. Rev. Dl, No.6, 1596 

Camenzind, M., 1990, in: Proc. of the IJfOth Symp. of I All, Heidelberg, 

June 19-23, 1989, Kluwer Academic Publishers, p. 413 
Celotti, A., Fabian, A.C., 1993, MNRAS 264, 228 
Cheng, K.S., Ho, C., Ruderman, M., 1986, ApJ 300, 500 
Clavel, J., et al., 1992, ApJ 393, 113 
Cohen, F., Piran, T., 1995, ApJ 444, L25 

Colgate, S.A., 1983, 18th Int. Cosmic Ray Conf., Vol.2, p. 230 
Daly, R.A., Guerra, F.J., Guijosa, A., 1996, to appear in: Energy Trans- 
port in Radio Galaxies and Quasars, eds. P. Hardee, A. Bridle, & A. 
Zensus (ASP Conference Series) 
Daugherty, J.K., Harding, A., 1982, ApJ 252, 337 

Dedenko, L.G., Zatsepin, G.T., Roganova, T.M., Fedorova, G.F., 1995 

Pis 'ma Zh. Eksp. Teor. Fiz. 61, No.4, 233 
Dermer, CD., Schhckeiser, R., 1993, ApJ 416, 458 
Dermer, CD., Miller, J. A., Li, M., 1995, ApJ 456, 106 
De Jager, O.C, Sleeker, F.W., Salamon, M.H., 1994, Nature 369, 294 
Dingus, B., et al., 1995, in :29*^ FSLAB Symposium "Towards the Source 

of Gamma- Ray Bursts", 25-27 April 1995, in press 
Drury, L O'C, 1983, Rep. Prog. Phys. 46, 973 
Dwek, F., Slavin, J., 1994, ApJ 436, 696 
Ferreira, J., Pelletier, G., 1993, A & A 276, 625 
Fichtel, C.F., et al., 1994, ApJS 90, 917 
Gallant, Y.A., Arons, J., 1994, ApJ 435, 230 
Gaisser, T.K., et al., 1993, Nucl. Phys. B 31, 399 

Gaisser, T.K., Halzen, F., Stanev, T., 1994, University of Madison Preprint 

MAD/PH/847 
Greisen, K., 1966, Phys. Rev. Lett. 16, 748 
Gould, R.J., Schreder, G., 1966, Phys. Rev. Lett. 16, 252 
Harris, D.F., CariUi, C.L., Perley, R.A., 1994, Nature 367, 713 
Hartmann, D.H., 1995, A & A Rev 6, 225 
Henri, G., Pelletier, G., Roland, J., 1993, ApJ 404, L41 
Hughes, P.A., AUer, H.D., AUer, M.F., 1989, ApJ 341, 54 



28 



Hurley, K., Dingus, B.L., Muherjee, R., et al., 1994, Nature 372, 652 
Hurley, K., Proc. of the Heidelberg workshop on TeV astrophysics, Eds. 

F. Aharonian & H. Volk, in: Space Science Reviews, Ed. De Jager, 

Kluwer, 1996, in press 
Jones, T.W., 1988, ApJ 322, 678 

Kamionkowski, M., 1992, in: Proc. of the Workshop on High Energy 

Neutrino Astrophysics, Honolulu, Hawaii, 23-26 March 1992, eds. V.J. 

Stenger et al.. World Scientific 
Kellermann, K.I., Pauliny-Toth, I.I.K., 1981, ARA&A 19, 373 
Kinney, A., 1994, in: Proc. of the First Stromlo Symposium: The Physics 

of Active Galaxies, eds. G.V. Bicknell, M.A. Dopita and P.J. Quinn, 

ASP Conf. Ser., Vol. 54, p. 61 
Lichti, G.G., et al., 1995, A & A, in press 

MacMinn, D., Primack, J., 1995, Proc. of the Heidelberg Workshop on 

TeV Astrophysics, Eds. H. Volk and E. Aharonian, Kluwer, in press 
Madau, P., Phinney, E.S., 1995, ApJ, in press 
Maisack, M., et al., 1995, A & A 298, 400 
Malkov, M., V61k, H.J., 1995, A & A 300, 605 
Mannheim, K., KriiUs, W., Biermann, P.L., 1991, A & A 251, 723 
Mannheim, K., Biermann, P.L., 1992, A & A 253, L21 
Mannheim, K., 1993a, A & A 269, 67 
Mannheim, K., 1993b, Phys. Rev. D 48, 5270 
Mannheim, K., Schhckeiser, R., 1994, A & A 286, 983 
Mannheim, K., 1995a, Astroparticle Physics 3, 295 
Mannheim, K., 1995b, A & A 297, 321 

Mannheim, K., Schulte, M., Rachen, J., 1995a, A & A 303, L41 
Mannheim, K., Westerhoff, S., Meyer, H., Eink, H.-H., 1995b, A & A, 
accepted 

Maraschi, L., et al., 1994, ApJ 435, L91 

Meisenheimer, K., et al., 1989, A & A 219, 63 

Meszaros, P., Rees, M.J., 1993a, ApJ 405, 278 

Meszaros, P., Laguna, P., Rees, M.J., 1993, ApJ 415, 181 

Meszaros, P., Rees, M.J., 1993b, ApJ 418, L59 

Meszaros, P., Rees, M.J., Papathanassiou, H., 1994, ApJ 432, 181 

Miyoshi, M., et al., 1995, Nature 373, 127 

Montigny, C. von, et al., 1995, ApJ 440, 525 

Parker, E.N., 1958, ApJ 128, 664 

Peebles, P.J.E., 1993, Principles of Physical Cosmology, (Princeton Univ. 
Press) 

Porcas, R.W., 1987, in: Proc. if a Workshop held at Big Bear Solar Ob- 
servatory, California, Oct. 28-30, 1986: Superluminal Radio Sources, 
eds. J. A. Zensus and T.J. Pearson, Cambridge University Press 

Pudritz, R.E., 1985, ApJ 293, 216 

Punch, C.W., et al., 1992, Nature 358, 477 



29 



Quian, S.J., et al., 1991, A & A 241, 15 

Quinn, J., et al., 1996, ApJ 456, L83 

Rachen, J. P., Biermann, P.L., 1993, A & A 272, 161 

Rawlings, S., Saunders, R., 1991, Nature 349, 138 

Roland et al., 1988, A & A 207, 16 

Rybicki, G.B., Lightman, A. P., 1979, Radiative Processes in Astrophysics, 

Cambrige University Press 
Urry, CM., Padovani, P., 1995, PASP 107, 803 
Sanders, D.B., et al., 1989, ApJ 347, 29 

Sigl, G., Schramm, D.N., Bhattacharjee, P., 1994, Astroparticle Physics 
2, 401 

Sikora, M., et al, 1987, ApJ 320, L81 

Sikora, M., Begelman, M.C., Rees, M.J., 1994, ApJ 421, 153 

Svensson, R., 1987, MNRAS 227, 403 

Stanev, T., et al., 1995, Phys. Rev. Lett. 75, 3056 

Sleeker, F.W., De Jager, O.C., Salamon, M.H., 1992, ApJ 390, L49 

Stenger, V.J., et al. (eds.), Proc. of the Workshop on High Energy 

Neutrino Astrophysics, Honolulu, Hawaii, 23-26 March 1992, World 

Scientific 

Tegmark, M., Hartmann, D.H., Briggs, M.S., Meegan, C.A., 1995, ApJ, 
submitted 

Thompson, D.J., et al., 1995, ApJS 101, 259 

van Putten, M.H.P.M., 1995, Nature, submitted 

V61k, H.J., Biermann, P.L., 1988, A & A 333, L65 

Walter, R., Fink, H.H., 1993, A & A 274, 105 

Webb, G.M., Axford, W.I., Terasawa, T., 1983, ApJ 270, 537 

Weekes, T.C., 1994, Bull. Am. Astr. Soc. 185, #60.01 

Zatsepin, G.T., Kuz'min, 1966, JETP Lett. 4, 53 



30 



